The energy scale behind the metalhc behaviors in low-density Si-MOSFETs 
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We show that the unexpected metalhc behavior (the so-called two-dimensional metal-insulator 
transition) observed in low-density Silicon metal-oxide-semiconductor field-effect transistors (Si- 
MOSFETs) is controlled by a unique characteristic energy scale, the polarization energy. On one 
hand, we perform Quantum Monte Carlo calculations of the energy needed to polarize the two 
dimensional electron gas at zero temperature, taking into account Coulomb interactions, valley 
degeneracy and electronic mobility (disorder). On the other hand, we identify the characteristic 
energy scale controlling the physics in eight different sets of experiments. We find that our ab-initio 
polarization energies (obtained without any adjustable parameters) are in perfect agreement with 
the observed characteristic energies for all available data, both for the magnetic field and temper- 
ature dependence of the resistivities. Our results put strong constraints on possible mechanisms 
responsible for the metallic behavior. In particular, there are strong indications that the system 
would eventually become insulating at low enough temperature. 



I. INTRODUCTION 

Until 1994, two dimensional systems were widely 
believed to be ultimately insulators^ as the presence 
of even a tiny disorder is enough to make the one- 
body electronic wavefunctions localized. In fact, un- 
til then, the experimental corpus was widely consis- 
tent with this paradigm^''^: the electrical resistance was 
found to increase as one lowered the temperature, in- 
dicating an eventual divergence at zero temperature. 
The observation of metallic behaviors in low-density 
Silicon metal-oxide-semiconductor field-effect transistors 
(Si-MOSFETs)^ hence came as a surprise and have been 
puzzling the community ever since^i^. The experiments 
were first acknowledged with skepticism but they were 
soon repeated and extended by several groups^r— , in- 
cluding in other materialsi^"— , so that it became clear 
that some intrinsic physics of the two-dimensional elec- 
tron gas was probed in those experiments. 

On the theory side, many different scenarios were pro- 
posed. In the weak disorder regime, perturbative calcu- 
lations performed in the eighties^Sr-.- showed the crucial 
role of spin and valley degeneracy and more recently, cal- 
culations in the limit of an infinite number of valleys sug- 
gested that electronic interactions could stabilize a non- 
Fermi liquid metallic phase^"^. Other theoretical works 
proposed that the metal-insulator transition could be due 
to an emulsion between a Wigner crystal and a Fermi 
liquid3^ or to some sort of a Wigner-Mott transitions^. 
It has also been argued that the metallic behaviors could 
be the signature of a superconducting phase^^"— . In con- 
trast to these rather " extreme" scenarios other works pro- 
posed a more conservative picture. This includes tem- 
perature dependent screeningSMS^ temperature depen- 
dent disorderi^iiSS or temperature dependent scattering 
between two different spin sub-bands split by a strong 
spin-orbit coupling'^^^'^^ . Percolation phenomenons were 
also claimed to play a role at the metal-insulator transi- 
tion itself'^^"^^ . The problem was also approached using 



numerical techniques and a delocalization effect by inter- 
actions was observed repeatedly^"—. 

The large spread in the above proposals, together with 
rather intense controversies between the conservative and 
the less conservative scenarios naturally lead to some 
confusion in the field so that no scientific consensus has 
emerged so far. The goal of this paper is twofold. First, 
we will perform a critical analysis of the existing experi- 
mental data and show that extremely strong constraints 
can be put on the possible theoretical scenarios. Second, 
following our recent proposal^, we will show that a sim- 
ple model with a minimum number of (well established) 
ingredients is enough to capture all the salient experi- 
mental facts. 

In section |TT1 wc summarize the main experimental 
facts with a focus on the characteristic energies that can 
be extracted from the data. In section Hill we introduce 
our minimum model to describe Si-MOSFETs, which is 
then used in section IIVI to calculate the polarization en- 
ergy Ep in presence of both disorder and electron-electron 
interactions. In section |Vl we come back to the exper- 
imental data and show that our calculated Ep agrees 
quantitatively (with no adjustable parameter) with the 
characteristic energies that can be extracted from either 
the magneto-resistance data or the resistance versus tem- 
perature data. Finally, we discuss the implications of the 
above findings for our own scenario^ (section |VI| as well 
as others (section IVlI|) . 



II. BASIC EXPERIMENTAL FACTS 

A. Main observations 

Let us quickly review the basic experimental facts that 
we want to understand. The first set of measurements 
that actually started the interest for these high mobil- 
ity MOSFETs is the behavior of the resistivity p{T) as 
a function of temperature (at temperatures lower than a 
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few Kelvin where phonon scattering no longer comes into 
play). The p{T) behavior is well understood in two limit- 
ing cases. At high density, the resistivity is rather small 
and depends only weakly on temperature, except for a 
weak negative dp/dT due to weak localization. At low 
density, the resistivity is much larger and one eventually 
reaches a clear insulating behavior where p depends very 
strongly (exponentially) on temperature {dp/dT < 0). 
Before the report made in Refi^, the common behavior 
observed in Si-MOSFETs, as well as in other two dimen- 
sional electron gases was a simple crossover between the 
high and the low density limit. Hence, the observation 
of an intermediate density regime with dp/dT > (here- 
after referred as property PI) came as a large surprise. 
In the left panels of Fig. [U we reproduced the data of 
three different experiments, including the original set of 
data from Kravchenko et al. One can observe a rather 
pronounced dp/dT > as the resistivity increases by 
a factor 10 between 100 mK and 4K. It was naturally 
surmised that in this regime, the resistivity could be ex- 
trapolated to a finite value at zero temperature, so that 
the system was in a metallic state, in contradiction with 
the prediction of scaling theory of localization!. 

The second important raw experimental feature is the 
dependence of the resistivity versus an in-plane magnetic 
field B: in presence of a magnetic field, the resistivity in- 
creases quickly before saturating (property P2). This is 
illustrated in the right panels of Fig. [1] for three different 
experiments. The increase of resistivity is comparable to 
the increase that one gets with increasing temperature, 
and dp/dT quickly becomes negative (insulating like). 

In our opinion, PI and P2 are the two main experi- 
mental features that require an explanation. There have 
been an important effort to establish scaling laws and 
critical exponents close to the metal-insulator transition, 
but we will not discuss this aspect. Indeed, once the 
metallic behavior PI is understood, there is no question 
that, as the system must be eventually insulating at low 
density, some sort of transition must occur. While this 
transition might be of interest, it is not fundamentally 
puzzling in itself. 
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FIG. 1: p as a function of T (left panels) and as a function of 
B (right panels), for different Ws, in various Si-MOSFETs ((a) 
from RelA, (b) from Ref", (c) from Ref.*^ (d) from Ref.^, 
(e) from Ref.— and (f) from Ref.^^). The red circles are the 
corresponding polarization temperatures Tp and polarization 
magnetic fields Bp obtained from the numerical simulations. 



B. First analysis 

One of the questions that many physicists had in mind 
when the first experiments'* came up was, what's new 
in those devices? Why have we not seen this behav- 
ior before? In fact, there was nothing qualitatively new 
about those MOSFETs. The real novelty lied in the fact 
that the samples were of extremely high mobility (for 
Si-MOSFETs) which allowed the experiments to be car- 
ried out at much lower densities than previously possible. 
Working at lower densities has three consequences: 

Property P3. The ratio of interaction over kinetic 
energy is large. This ratio — ra* / {^-Kt^tt? yjixus) 
(where ris is the electronic density, m* = 0.19 rrie the ef- 
fective mass, e the electron charge, eg the dielectric con- 



stant and e = 7.7 the relative dielectric constant) takes 
values between 5—10 for the metallic region, hence the 
system is intrinsically correlated. 

Property P4. The effective disorder seen by the elec- 
trons, as parametrized by l/kpl {kp Fermi momentum 
and I mean free path) is paradoxically not smaller: in- 
deed the high mobility of the samples is compensated by 
the fact that the density is lower and therefore kpl oc Us 
remains of the order of a few units. Note that contrary to 
ris (hence Ts) which can be almost directly measured with 
Shubnikov-de-Haas measurements, there is no generic 
way to measure the disorder. At high density, r-g <C 1 
and the interaction effects can be neglected, so that kpl 
is simply related to the conductance g of the system, 
g — (2e^ /K)kpl. In this limit, we estimate the product 
r\ = rs\fkpl which should remain constant as one lowers 
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the density (as oc l/y^nj and kpl (x Us). Hence, rj is 
estimated at large density by the non-interacting formula 
77 = y/JIe'^/^m* /{4:TThp/^eoe). For real samples where dis- 
order have intrinsic characteristic lengths, the mobility 
fi may have a density dependence. Typical values of rj 
are of the order of 20. In Fig. [21 we have collected the 
trajectories in the versus l/y/kpl phase diagram for 
ten different samples: as one lowers the density, the sam- 
ples follow a straight line l/y/kpl = rg/rj. The regions 
where the samples have a metallic behavior are plotted 
in blue while the insulating behaviors appear in dashed 
red. When comparing the samples where the metallic 
behavior is observed with the older samples that did not 
show such a behavior (the two on the left), one indeed 
observes that the metallic behavior is observed at larger 
Ts (P3) but in the same regime of disorder (P4). 




FIG. 2: Experimental rs versus I/V^f' phase diagram of the 
Si-MOSFET, constructed from ten samples studied in the lit- 
erature. Each sample corresponds to a straight line in the 
diagram (dotted lines). Blue (dashed red) refers to the den- 
sity range where metallic dp/dT > (insulating dp/dT < 0) 
behaviors are observed. This diagram should be compared 
with Fig.4 of Ref.*^ 

Property P5. The Fermi energy Ep decreases so 
that the ratio T/Tp increases (where Ep — kTp — 
Trh^ns/2m* , Tp is the Fermi temperature). Typical ex- 
perimental values are quite large, T /Tp w 20% around 
the inflexion point of the p(T) curves as can be seen 
in Fig. [l](c). This last point, raised in Ref.-^ has been 
overlooked by many authors and is however crucial: the 
metallic behavior happens for temperatures which are of 
the same order of magnitude as the Fermi temperature. 
Hence, the usual "low energy" paradigm where one only 
looks at electrons close to the Fermi surface does not hold 



in those devices. 

It is actually the purpose of this paper to identify the 
characteristic energy that controls the behavior of p{T) 
and p{B). As we shall see, both are controlled by a 
unique energy scale, the polarization energy Ep of the 
Si-MOSFET. 



III. A MINIMUM MODEL FOR SILICON 
MOSFETS 

We now introduce the model that we will use for 
our calculations. It includes Coulomb repulsion, disor- 
der, spin and a twofold valley degeneracy (present in Si- 
MOSFETs but not in their GaAs counterparts, Silicon 
being an indirect gap semi-conductor). This model is 
minimum in the sense that the first three ingredients are 
already known to be relevant experimentally, as discussed 
in the above section, while the fourth one is known to be 
present and we will show that it plays an important role 
for quantitative predictions. The system contains iV ^ 1 
particles in a lattice made of x Ly sites, the particles 
being equally split up into the two degenerate valleys, 
with spin up and iV^ spin down. The spin configura- 
tion defines the system polarization 
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The Hamiltonian for the 4-components plasma reads. 



H 



{fr')(7 



ra^r' a' 



(2) 

where ct^ and cpa are the usual creation and annihilation 
operators of one electron on site r with inner degree of 
freedom cr, the sum X](fr')(T restricted to nearest neigh- 
bors and np„ = ct cpa is the density operator. The inter- 
nal degree of freedom a = 1 ... 4 corresponds to the spin 
and valley degeneracy degrees of freedom. The disorder 
potential is uniformly distributed inside [— VF/2, W^/2]. 
t is the hopping parameter and U is the effective strength 
of the two body interaction Vp. To reduce finite size ef- 
fects, is obtained^ from the bare Coulomb interaction 
using the Ewald summation technique. 
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L 



■Erfc(fcc|f +L|) 



(3) 



+ 7^ E ^Erfc(|ie|/(2fc,)) cos(7? • f) . 
\K\ 

In the previous equation, k^ is a (irrelevant) cut off. 
The vector L takes discrete values L = {nxLx,nyLy) 
with Ux and Uy integer numbers. The vector K also 



takes discrete values, K = {- 



and {nx,ny) ^ 
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(0,0). The complementary error function is defined as 
Erfc(a;) = ^ e~* dt. We work at small filling fac- 
tor u = N/{LxLy) <^ 1, where we recover the continuum 
limit. The two dimensionless parameters that control the 
physics read, 



(4) 



Ts = U/{2ty/^), kpl = iSiTiyr/W^ 



IV. QUANTUM MONTE-CARLO 
CALCULATION OF THE POLARIZATION 
ENERGY. 

This section is devoted to the numerical calculation of 
the polarization energy of the 4-components plasma, as 
a function of and kpl. We compute the ground state 
energy per particle E(j>) of the system with the Green 
function Quantum Monte Carlo (QMC) technique in the 
Fixed Node approximation (see Refj^ for the numerical 
method). Particular care was given to the extrapola- 
tion to the thermodynamics limit {N ^ 1) as well as 
the continuum limit (the algorithm being implemented 
on a lattice). All the fits given below, as well as data 
points, are given with a precision better than ±0.02Ep 
(for kpl > 1.5) and ±OMEp (for 0.3 < kpl < 1.5). The 
data have been averaged on 50 to 200 samples for the 
strongest disorders. 

At very high density, i.e. no disorder {1/kpl = 0) and 
interaction (r^ = 0), it is straightforward to show that 
the energy depends on the polarization p in a quadratic 
manner: E{p) = E p{l + p^) / 2 . In fact, we find that such 
a quadratic dependence is verified with good precision 
for a wide range of parameters < r<, < 10 and kpl > 1, 
corresponding to the disorder regime of the experiments, 
as shown in the upper left panel of Fig. |31 We found 
(small) deviations from the quadratic behavior for rather 
large disorder in the localized regime (see triangles at 
1/kpl = 2.65 and = 6). Hence we find with good 
precision that. 



E{p) - E{0) + Ey , 



(5) 



where the coefficient Ep{rs, kpl) is the polarization en- 
ergy of the system and is proportional to the Fermi en- 
ergy. Without electron-electron interaction and disorder 
Ep — Ep/2, and as interaction is turned on, Ep decreases 
(indeed a ferromagnetic state tends to minimize Coulomb 
energy as antisymmetric orbital wave-functions have a 
lower probability for two electrons to be close to each 
other than symmetric ones). 

In presence of an in-plane (i.e. purely Zeeman) mag- 
netic field B, Eq.® turns into £'(p) = £^(0) + Epp^ - 
g^BBp/2 [g — 2 Lande factor and /zb Bohr magneton). 
Introducing Bp = AEp/{gfiB), the minimum of energy is 
found for a polarization p* given by. 
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P 



B/Bp 
1 



B <Bp 
B> B„ 



(6) 
(7) 
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2 


3 


4 


5 


A, 


27.93 


9.83 


56.5 


46 


1.77 


0.019 



TABLE I: Parameters Ai of Eq. (O and Q. 



Hence, the spin susceptibility x — nsl2{dpldB)B=ti is 
directly related to the field Bp at which the polariza- 
tion saturates. Defining the non-interacting susceptibil- 
ity Xo = g^J.Bm* /{2Trh'^) , one finds Ep/Ep = 2x/xo- In 
addition to the characteristic magnetic field Bp, Ep is also 
associated with the characteristic temperature kTp = Ep 
at which the polarized excited states get significantly 
populated. 

In the lower panels of Fig. [31 we plotted our Ep data 
as a function of disorder l/y/kpl (left) and interaction 
Ts (right). Without disorder (squares in the lower right 
panel), our results are in close agreement with previous 
Diffusion Monte Carlo calculations performed in Refi^ 
(dashed line). As expected, the polarization energy Ep 
decreases with interaction r^. However, contrary to the 
single- valley case where Ep becomes extremely small and 
could even become negative (ferromagnetic instability), 
we find that Ep is always positive up to Ts = 20 where 
Ep = O.lSEp. 



We find that disorder tends to increase Ep for very 
small interactions Vg < 0.25. At — 0, the disorder 
correction to Ep is very well described by the second or- 
der perturbative correction Ep/Ep = 1/2 + \og2/ {■Kkpl) 
(dashed line in lower left panel of Fig. However, as 
soon as small interactions are switched on (rj > 0.3), dis- 
order tends to decrease Ep. We also find that the second 
order perturbative correction oc 1/kpl is valid only for 
very tiny disorder 1/kpl < 0.04 above which the correc- 
tion becomes proportional to 1/y/kpl. In order to make 
our set of data easily retrievable, we use the following fit, 
vafid for our entire data set for 0.25 < r., < 10, 



Ep{kpl,rs) = E;\rg) + 



Ep+A^Ep, (8) 



where Ep^ is the polarization energy of the clean system 
fitted (for < < 20) with the Fade formula: 



Aq + Ai Ts 



A2 + A:irg + Aa r1 



Ep 



(9) 



The values of fitting parameters Ai are listed in Table HIl 
The parameter a is plotted in the upper right part of 
Fig. |3] and is equal with good precision for rs > 2 to 
a = —0.1. The small parameter A^ roughly accounts for 
the fact that at very low disorder, the disorder-correction 
to Ep/Ep is not linear but quadratic in l/^/kpl. We 
note that Eq.® and ^ are merely a compact way of 
describing our numerical data. 
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FIG. 3: Polarization energy. Upper left panel: E{p) as a 
function of p'^, for l/kpl — 0, fs = (diamonds), l/kpl = 0, 
Vs — 10 (circles), l/kpl = 0.66, = 2 (squares) and l/kpl = 
2.65, rs — 6 (triangles). Straight lines are linear fits. Upper 
right panel: a as defined in Eq.(j8} as a function of rs- Lower 
left panel: Ep as a function of disorder l/^/kpl, for Ts = Q 
(circles), ra = 0.25 (left triangles), Vs = 0.5 (down triangles), 
Ts = 1 (diamonds), rs = 2 (up triangles), rs = 5 (squares) 
and rs = 8 (right triangles). Dashed line stands for Ep/Ep ~ 
1/2 + log2/(7rfcF0- Lower right panel: Ep as a function of rs 
for l/kpl ~ (squares), l/kpl = 0.66 (circles), l/kpl ~ 1.49 
(diamonds) and l/kpl = 3.51 (triangles). Dashed line is the 
result of Ref.^^ without disorder. In both lower panels, the 
lines are fits given by Eq.® for rs > 0.25. All energies are in 
unit oi Ep- 

V. COMPARISON OF THE POLARIZATION 
ENERGY Ep WITH THE EXPERIMENTAL 
CHARACTERISTIC ENERGY SCALE. 

Fig. [T] shows raw experimental data p{B) and p{T) on 
low density Si-MOSFETs from six different experiments, 
including the original set of data from Kravchenko et 
al. For each sample, we estimate the parameter 77 from 
the peak mobility ^p"^'^ obtained at high density (hence 
probably overestimating rf) which sets how y/kpl — rj/rs 
evolves with density. Then, for each density, we calculate 
the polarization energy Ep according to Eq. ([8]) and ([9|) . 
Disorder corrections to Ep are rather small for those high 
mobility samples, so that our results are dominated by 
Ep{rs) and not very sensitive to our estimate of 77. For 
each curve p{B) {p{T)), we draw a circle at the calculated 
saturation field Bp (characteristic temperature Tp). The 
calculations are performed without any adjustable pa- 
rameters and we find an extremely good match between 



Bp (Tp) and the field (temperature) at which the resis- 
tivity saturates (has its inflexion point). 



A. Magneto-resistance experiments 

Describing the full p(i3)-dependence is a complicated 
task. However, as the only effect of an in-plane magnetic 
field is to polarize the sample through Zeemann coupling 
(the orbital effect of an in-plane field can be neglected 
with good approximation), one expects that when the 
polarization saturates, the resistivity also saturates^, so 
that p{B) provides a direct measurement of Bp. 

The experimental data of resistivity as a function of 
in-plane magnetic field B (at low temperature) shown 
in panels (d), (e) and (f) of Fig. [T] correspond to 
three samples with mobility ^p^^"*-' = 25000 cm^/Vs^, 
41000 cmVVs^ and 20000 cmVVs^. We find that the 
magnetic field B^at at which the resistivity saturates is in 
extremely good agreement with our calculated Bp. This 
agreement between the experimental B^at and the calcu- 
lated Bp can be considered as a validation of the mini- 
mum model that we have used as well as of the accuracy 
of the QMC method. 

Let us now discuss a publication^^ which claimed to 
observe a divergence of the spin susceptibility x- Such 
a finding would contradict our calculation as x cx 1/Bp, 
and although Bp decreases as one lowers the density, we 
did not observe any divergence in the numerical calcu- 
lations (one should however keep in mind that in the 
lower right panel of Fig|31 the energies are in unit of 
Ep which itself is proportional to the electronic den- 
sity Us). The corresponding data are shown in the 
right panel of Fig. |4] (our Bp-, full line, experimental 
Bsat- dashed line). The experimental B^at were obtained 
in Ref.i^' from rescaling the in-plane magneto-resistivity 
data inP""'' = 30000 cm^/Vs) at smah field. Here, we find 
a good, but not very good, agreement between Bp and 
Bsat- To understand the discrepancy, let us come back 
to the raw p{B) experimental data for two values of the 
density corresponding to the blue squares (experiments) 
and red circles (our numerics) . Those raw data are shown 
on the left panel of Fig. |4] where we also reproduce our 
predictions as well as the value of Bsat obtained from 
the rescaling procedure. Clearly, the rescaling procedure 
underestimates Bgat as the resistivity has not saturated 
yet at this value while our calculated Bp matches pre- 
cisely the end of the resistivity rise. We find that the 
agreement between the experiments and the numerics is 
actually extremely good and that the rescaling procedure 
only give approximate estimates of Bsat- As the diver- 
gence of X claimed in Ref.^^ is extremely sensitive to this 
procedure, we conclude that there is probably no diver- 
gence of the susceptibility as defended in Refsj ^^'^^ but in 
contradiction with Ref.— . New thermodynamic measure- 
ments are in qualitative^ and quantitativeSS agreement 
with our numerics (hence we disagree with the conclu- 
sions of Ref.^° which claim that there is a divergence of 
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B (Tesla) 

FIG. 4: Left panel: p as a function of B at different Us, from 
Ref.— . Right panel: experimental Beat from Ref.— (dashed 
line), QMC Bp (full line) and experimental Bp from Refi^ 
(green squares) as a function of r^. In both panels, red points 
correspond to QMC Bp and blue squares to Bsat from Ref.^^, 
at two fixed densities. 

the spin susceptibility). 

Finally, we add a last data set that corresponds to the 
measures made in Ref;^ (green squares in right panel 
of Fig. U) from the analysis of Shubnikov de Haas oscil- 
lations, at a similar mobility [pP^°-^ = 24000 cm^/Vs). 
We find that they agree perfectly with our calculations. 
Overall, we have a quantitative agreement between Bp 
and Bsat for five different experiments indicating that 
our minimum model and QMC calculations adequately 
describe the experimental situation. 



B. Resistance versus temperature experiments 

We now turn to the original puzzle, the resistivity ver- 
sus temperature p{T) experiments. A first remark that 
was originally made in Refi^ is that the large increase 
of resistance takes place at large temperatures compare 
to usual transport experiments. This can be seen for in- 
stance in Fig. [T](c) where temperature has been rescaled 
with respect to the Fermi temperature Tp . One finds that 
the infiexion point of p{T) takes place for T/Tp = 0.2. 
Refs.'^^'- ^ attributed this large energy scale to some semi- 
classical effect. Here however, we argue that this char- 
acteristic scale is the polarization temperature Tp. On 
a qualitative level, it is not surprising: this energy scale 
is equal to a fraction of the Fermi temperature (0.5 Tp 
without interaction and disorder and less in their pres- 
ence) and we have already showed that it controls the 
magneto-resistance behavior. 

More quantitatively, in panels (a), (b) and (c) of Fig. [1] 
we plot various experimental data of resistivity p as a 
function of temperature T, at zero magnetic field and 
at various electronic densities, for three Si-MOSFETs 
with respectively fj.P""'' = 30000 cm^Vs (Ref.-^), 
41000 cmVVs (Ref^) and 19600 cmVVs (Ref4^). We 
have superimposed the circles corresponding to Quantum 
Monte Carlo polarization temperatures Tp. The typical 



temperature where one observes the change of resistivity 
associated with the metallic behavior matches our calcu- 
lated Tp. This is a strong indication that the polarization 
energy is the unique energy scale that controls the physics 
of the metallic behavior. 



VI. A MINIMUM SCENARIO FOR THE 2D 
METAL-INSULATOR TRANSITION 

We now have a good understanding of what happens 
on the horizontal of the p{T) and p{B) curves. To gain a 
full understanding of properties PI and P2, we need to 
understand what happens on the vertical axis, i.e. why 
does the resistivity increase and not decrease upon in- 
creasing temperature or magnetic field. In other words, 
in order to account for both PI and P2, it is necessary 
and sufficient to show that polarized states have a higher 
resistivity than non polarized one (hereafter referred as 
property P6). Indeed, it is necessary, since the p{B) 
curves show that the polarized system (high field) has a 
higher resistivity than the non polarized one (zero field). 
It is also sufficient as, once P6 is established, it natu- 
rally follows that when one will increase temperature in 
a range around Tp, the highly resistive excited states will 
be significantly populated and the overall resistivity will 
increase. 

In a previous publication^, we have performed a sys- 
tematic study of the interplay between Anderson localiza- 
tion and electron-electron interactions in Si-MOSFETs. 
We found that upon increasing interactions, the localiza- 
tion length of the non polarized ground state strongly in- 
creases (in absolute value, but in particular with respect 
to the polarized excited states), so that one naturally ac- 
counts for P6. This result is also shown in the right panel 
of Fig. [5] where we compare as a function of the local- 
ization lengths ^ (in unit of the average distance between 
electrons a) computed at zero temperature for the non- 
polarized electron gas and for the fully polarized one. At 
very weak interaction (t's < 0.5), the polarized system is 
less localized than the non-polarized one. Indeed, polar- 
izing the system raises the Fermi level Ep (by a factor 
2), hence also kpl, hence ^ which depends exponentially 
on the latter. On the contrary, at stronger interaction 
(^s }t 0.5), the situation is reversed and the polarized 
system is more localized than the non-polarized one. We 
checked that we also get such an inversion of the localiza- 
tion lengths between the non-polarized ground state and 
partially polarized states, so that we indeed recover P6. 
In Fig. 4 of Ref.^^ we presented the effective phase dia- 
gram constructed out of this mechanism which should be 
compared with the experimental one presented in Fig. [51 
The agreement is semi-quantitative, meaning that our 
scenario also captures properties P3 and P4 (as well as 
P5 discussed above). Hence, we find that the minimum 
model discussed in this paper is enough to account for 
all the relevant experimental facts discussed in the be- 
ginning of this paper. 
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FIG. 5: Left panel: a/^ as a function of Xj^kpl for the non- 
polarized Si-MOSFET. The symbols correspond to different 
strengths of interaction: from left to right, rs = (circles), 
Ts = 0.1 (left triangles), = 0.25 (full squares), — 0.5 
(empty diamonds), rs = 1 (up triangles), then from right to 
left, Ts = 2 (down triangles), rs = 4 (right triangles), rs = 6 
(empty squares) and Ts — 10 (full diamonds). Right panel: 
a/5 as a function of rs for the non-polarized (NP, full lines) 
and polarized (P, dashed lines) Si-MOSFET, at kpl = 1.51 
(two lower green curves) and kpl = 0.89 (two upper blue 
curves). The arrows indicate the interaction strength rs at 
which the polarized excited state becomes more localized than 
the non polarized ground state. In both panels, the lines are 
given by Eq. fTD]), fTT]) and ([T^. 

We end this section with a discussion of the existence 
of "true" metal at zero temperature, i.e. does the local- 
ization length diverge in the region where the metallic 
behavior is observed. We emphasize that in light of the 
above discussion, this is a rather academic question as 
what has actually been measured experimentally corre- 
sponds to rather high energy physics. In Ref.*^, we have 
claimed that there probably does not exist a " true" metal 
at zero temperature. Our main argument was the ab- 
sence of any visible deviation to the one-parameter scal- 
ing theory^ in presence of interactions, despite a huge 
delocalization effect. Here we give another argument. 
In the left panel of Fig. [5l we plot a/^ as a function 
of the disorder strength 1 / ^/kpl for various values of rg . 
Without interaction^, the localization length of the two- 
dimensional electron gas depends exponentially on kpl: 
£,/a = {kFl/B)exp{Ckpl). We find that in the pres- 
ence of interactions our entire set of data can be fitted 
by the same law upon renormalizing the B and C pa- 
rameters. Hence, while the localization length can be- 
come large in the presence of interaction, it always re- 
mains finite and the system remains ultimately an insu- 
lator. In practice this insulating behavior should appear 
at extremely low temperature when the phase coherence 
length becomes larger than the localization length. Such 
a metal-to-insulator re-entrance as a function of temper- 
ature has been observed at high densitj^i^ and at low 
density, in the vicinity of the so-called metal-insulator 
transition, in a GaAs two-dimensional hole gas^ and in 
Si-MOSFETs?fi. 

Below we provide the fitting parameters that account 



for our localization length data. We have 

^{kpl, r,) = exp {-C{r,)kFl) , (10) 

1^ kpl 

where two different sets of parameters are needed for the 
non-polarized {B^^ ,C^^) and the polarized {B^,C^) 
case. In turn, the parameters B and C are well fitted by: 

B{rs) = bo + bly/r^ + b2rs + b3rs^/r^ 

+ 64r2 + 55r2v?r+66'^' (11) 

and: 

C{rs) = cq + ciy/r^ + C2rs + csrs^/r^ 

+ C4r^s+C5r^sV^s+cerl, (12) 

where the values of the bi and Ci fitting parameters are 
given in Table IhI for the non-polarized system (6f ^ and 
cf'^) and for the polarized one {bf and cf ), with the 
convention that kpl always refers to the Fermi level of the 
non-polarized gas. Eq. pUj) supports the idea that the 
two-dimensional electron gas behaves as a Fermi liquid 
(at least for moderate interactions < 10) and that the 
effect of electron-electron interactions is to renormalize 
its effective characteristics. 



VII. CONCLUSION 

Let us discuss a few theoretical scenarios that have 
been proposed by other authors. In light of the findings 
of this paper, we believe that models whose characteristic 
energies do not involve the polarization energy should not 
be applied to those experiments. This includes for exam- 
ple the percolation scenario proposed in Ref.^^ (although 
percolation probably plays an important role close to the 
insulating region in GaAs heterostructures) . This also 
includes the proposal related to Wigner crystal Refi^ 
(which we dismissed in Ref.^'^ on other grounds). Refi^ 
has some overlap with our scenario (on the role of polar- 
ization energy in particular) but relies on the existence of 
some form of local order (which we did not observe in our 
simulations) and does not include the presence of disor- 
der. Extrinsic models as the temperature dependent dis- 
order scenario proposed in Refsi^i^^ can also be excluded 
on the same ground (but these authors were the first to 
recognize the presence of a high energy scale in the exper- 
iments). An important theoretical advance was made in 
Refi^ where the authors studied the interplay between 
interaction and disorder in the diffusive limit, in the limit 
of an infinitely large number of valleys. We believe that 
those authors correctly pointed out the important role 
of valley degeneracy in the experiments, at the origin of 
the large differences between the behaviors observed in 
Silicon MOSFETs and in other heterostructures. How- 
ever, the limit of an infinite number of valleys artificially 
increases the role of electron-electron interactions. In 
particular those authors find that even a very weak in- 
teraction will drive the system toward a non Fermi-liquid 
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i 





1 


2 


3 


4 


5 


6 




0.191 


-0.325 


0.24 


-0.0526 


0.0035 










0.136 


-0.09093 


-0.03197 


0.18232 


-0.17749 


0.07549 


-0.010827 


NP 


0.84 


-0.177 


1.094 


2.0616 


-1.6784 


0.3027 





cf 


1.68 


-0.712 


2.418 


-1.1874 


0.1594 









TABLE II: Parameters h and d of Eq. ^ and ^ for the non-polarized Si-MOSFET (fcf ^ and cf ^) and for the polarized 
one (6f and cf ). 



fixed point. Those results, which imply a complete break- 
down of one parameter scaling theory are in contradic- 
tion with our numerical results^. More importantly, the 
corresponding analysis of the experiments^!^ is in con- 
tradiction with the present analysis of the characteristic 
energy scales involved in these systems. 

In summary, we have considered a minimum model 
of Si-MOSFETs taking into account Coulomb repulsion, 
spin, valley degeneracy, and disorder in a non perturba- 
tive way. Our chief result is that the polarization energy 
Ep calculated for this model is in quantitative agreement 
with the characteristic energy scale that controls the 
metallic behavior of those high mobility Si-MOSFETs. 
Beside spin. Coulomb repulsion plays a crucial role here 
as it decreases Ep by a factor 2 — 3 (see FiglS]). The 
presence of valley degeneracy is also very important as 
it also decreases Ep by a factor 2. On the other hand, 
disorder only gives corrections of 10 — 30% to Ep for 
the samples considered in this survey. We note that the 
polarization temperature is the "crossover temperature" 
that was forseen in Ref.^^. However, we find that this 
"crossover temperature" is an intrinsic property of the 
electron gas, and does not involve any additional extrin- 
sic ingredient. A simple corollary of our chief result is 
that the non-polarized system has to be a much better 
conductor than the polarized one (to explain the sign of 
dp/dT and dp/dB in the metallic region). We have ar- 
gued before^ that our minimum model accounts for this 
counter- intuitive point (P6), as the interplay between 
Coulomb repulsion and Anderson localization depends 
strongly on the polarization. These findings imply that 



our "necessary" model (in the sense that all its ingredi- 
ent are known to be present and a priori relevant) is also 
"sufficient" to capture the essential features of the metal- 
lic behaviors. Another consequence of our model is the 
apparent failure of one parameter scaling. For a given 
polarization^ we could not find any deviation to the (one 
parameter) scaling theory of localization^. However, at 
finite temperature excited states of different localization 
lengths come into play so that the physics is no longer 
controlled by a single parameter. In this work, we have 
focused on Si-MOSFETs, but at the qualitative level, 
many statements also apply to other materials like GaAs 
heterostructures. In our view, the situation in the latter 
is made a bit different by the conjunction of three ele- 
ments. First the absence of valley degeneracy makes the 
delocalization effect of the ground state much less effec- 
tive. Second, the present of the doping layer close to the 
gas induces some long range disorder which can mask the 
local physics. Last, extremely high mobility are available 
for these systems, making it possible to study almost bal- 
listic samples. 
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